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O ' Abstract 

(N : 

In this paper, we consider an initial boundary value problem for Maxwell's equa- 
P_( ! tions. For this hyperbolic type problem, we derive guaranteed and computable upper 

' bounds for the difference between the exact solution and any pair of vector fields in 

the space-time cylinder that belongs to the corresponding admissible energy class. 
For this purpose, we use a method suggested in |20] for the wave equation. 
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1 Introduction 

In this paper, we derive computable upper bounds for the distance between the exact 
solution {E, H) of an initial boundary value problem for Maxwell's equations, which 
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have in their second order form a hyperbohc nature, and any pair of vector fields {E, H) 
belonging to the admissible energy class of the considered problem. As our techniques rely 
on second order methods and the Maxwell system decouples in its second order version 
for the electric field E and the magnetic field we focus on E in our analysis. The 
vector field E can be considered as an approximation of E computed with the help of 
a numerical method. In other words, we deduce nonnegative functionals B (also called 
error majorants or upper bounds) that depend only on E and known data (coefficients, 
domain, right hand side and boundary data) and satisfy the following properties: 

1. £{E -E)< B{E) for all admissible E. 

2. B{E) = OiSE = E. 

3. B{E) if S{E-E) 0. 

Here, £^ is a suitable error measure of the energy of the system defined on a space-time 
cylinder (e.g., a L^-energy norm). 

Such Functionals B provide an explicit verification of the accuracy of approximations. 
Indeed, we see that B{E) is small then E belongs to a certain neighborhood of the 
exact solution. Moreover, B vanishes only at the exact solution E. The third property 
shows that the majorant B possesses the continuity property with respect to all sequences 
converging in the topology induced by the energy norm S. 

Estimates of such a type (often called functional a posteriori estimates) can be derived 
by at least two methods. The first method is based on variational techniques and appli- 
cable for problems that admit a variational statement. By this method a posteriori error 
estimates were derived in [HI [T5] and many other publications (see [5] for a systematic 
overview). Another method is based upon the analysis of the integral identity (varia- 
tional formulation) that defines the corresponding generalized solution. This method was 
suggested in [16], where it was also shown that for linear elliptic equations both meth- 
ods (variational and nonvariational) lead to the same estimates. Later the nonvariational 
method was also applied to nonlinear elliptic problems and to certain classes of nonlin- 
ear problems in continuum mechanics (e.g., for variational inequalities [H [21 US]) and to 
initial boundary value problems associated with parabolic type equations [T7]. A conse- 
quent exposition of the 'nonvariational' a posteriori error estimation method is presented 
in the book [19]. Analogous estimates have been derived for elliptic problems in exterior 
domains as well [TO] . 

In this paper, we are concerned with an initial boundary value problem for Maxwell's 
equations. For the stationary version of this problem, functional a posteriori estimates 
have been derived earlier in [11] for bounded domains. However, the hyperbolic Maxwell 
problem essentially differs from the stationary case and the estimates are derived by a 
new technique. The derivation method is also based on the analysis of a basic integral 
relation but uses a rather different modus operandi. The reason for this lies in the specific 
properties of the respective differential operator involving second order time and spatial 
derivatives with opposite signs. We overcome the difficulties arising due to this fact with 
the help of a method suggested in [20] for the wave equation, which is closely related. 
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and deduce computable upper bounds for the distance to the exact solution measured in 
a canonical L^-energy norm. 

Our main results are presented in Section |3] by Theorems [5] and [8|, which provide 
computable and guaranteed majorants for the error measures (13. 2p and (13. 3 p of the electric 
field E. These first (and simplest) majorants are derived under stronger assumptions on 
the approximation. They can be used if the approximation E possesses extra regularity, 
which sometimes may be difficult to guarantee in many numerical schemes. In Section H] 
we prove corresponding results under weaker assumption on the approximation E, which 
are free of these drawbacks, but have a more complicated structure. Finally, in section [5] 
we estimate the error for the approximation of the magnetic field H as well and thus the 
error of the approximation of the full solution {E, H). 

We note that the respective functionals generate new variational problems, where exact 
lower bounds vanish and are attained only on the exact solution. In applied analysis, the 
functionals can be used for a posteriori control of errors of approximate solutions obtained 
by various numerical methods. 



2 Basic problem 

Let f2 be a domairQ in M.^ with Lipschitz continuous boundary T := dD. and correspond- 
ing outward unit normal vector by u. Furthermore, let T > 0, I := {0,T) as well as 
Qf := {0,t) X Q and Fj := (0, t) x F for alH > the space-time cylinder and cylinder bar- 
rel, respectively. We consider the classical initial boundary value problem for Maxwell's 
equation: Find vector fields E and H (electric and magnetic field), such that 

dtE - £-^cnr\H = F in fiy, (2.1) 

dtH + fi'^cmlE = G in (^r, (2.2) 

ux E\r = on Fy, (2.3) 

^(0) = ^(0, ■ ) = ^0 infi, (2.4) 

H{0) = H{0, ■) = Ho in (2.5) 

Here e and /i denote time-independent, real, symmetric and positive definite matrices with 
measurable, bounded coefficients that describe properties of the media (dielectricity and 
permeability, respectively). For the sake of brevity, matrices (matrix- valued functions) 
with such properties are called 'admissible'. We note that the corresponding inverse 
matrices are admissible as well. 

Remark 1 The underlying domain Q may be bounded or unbounded. Contrary to the 
stationary cases, i.e., static or time-harmonic equations, the Sobolev spaces used for the 
solution theory of the Cauchy problem do not differ whether the domain is bounded or not. 
For instance, in exterior domains one has to work with polynomially weighted Sobolev 
spaces what naturally would lead to weighted error estimates as well. See |3 0, 
for a detailed description. 



*, i.e., a connected open set, 
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By L^(f2) we denote the usual scalar L^-Hilbert space of square integrable functions 
on Q and by H™(i7), m G N, the usual Sobolev spaces. H(f2) denotes the Hilbert space 
of real-valued L^-vector fields, i.e., L^(f2,M^). For the sake of simplicity, we restrict our 
analysis to the case of real-valued functions and vector fields. The generalization to 
complex-valued spaces is straight forward. Moreover, we define 



H(curl,fi) := |$ G H(fi) | curl$ G H(l])} , H(curr,fi) := C°°(fi), 

where the closure is taken in the natural norm of H(curl, Q). The homogeneous tangential 
boundary condition (12. 3p is generalized in H(curr,fi) by GauB' theorem. Equipped with 
their natural scalar products all these spaces are Hilbert spaces. 

To formulate and obtain a proper Hilbert space solution theory for the latter Cauchy 
problem, we need some more suitable Hilbert spaces. We set 

n{Q) := H(fi) X H(fi) 

as a set and equip this space with the weighted scalar product 

H(n)xH(Q) H(n) ' 



where 



A : = 



e 
/i 



For the sake of a short notation, we will write for domains S C 

II ■ IIh •= II ■ IIl2(h,m<') ' (■,■)=•=(■'■ )l2(H,IR'^) 



and for suitable matrices A 

Furthermore, we introduce the linear operator 

Ma : V{Ma) C H{n) H{n), ($, ^) ^ iA-^M{^, ^) 

putting 

V{Ma) := H(curr,fi) x H(curl,n), M : 




curl 



curl 




Then, a solution of the Cauchy problem fl2.1l) - fl2.5p is to be understood as a solution of 
the Cauchy problem 



idt-iMA)iE,H) = {F,G), 
iE,Hm = iEo,Ho) 



(2.6) 
(2.7) 
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Utilizing a slight and obvious modification (variation of constant formula) of [U The- 
orem 8.5], the Cauchy problem fl2.6p -( ET|) has unique solution for all T (we may also 
replace the interval / by M) by spectral theory since M\ is self-adjoint. The spectral 
theorem suggests 

{E,H){t) = exp{itMA){Eo,Ho)+ [ exp{i{t - s)Ma){F,G){s) ds, t el 



as solution. We get: 

Theorem 2 Let (F, G) G l\l,H{n)) and {Eq, Hq) G H{n). Then, the Cauchy problem 
fl2.6p - fl2.7l) is uniquely solvable in 

(i) C^IM^)); 

(ii) C°(l,V{MA))nC(l,n{n)), if additionally 

{F,G) G 1\I,V{Ma)) n C^(l,nm and {E,,H^) G P(Ma); 

(iii) C^{I,V{Ml))r}C{I,V{MK))r]C(!,H{n)), if additionally 

{F,G) G \}{I,V{Ml))nQ.\l,V{MA))r\Q\l,U{n)) and {Eo,H,) G V{Ml). 

Here, {E, H) G V{Ml), if and only if 

{E,H),{e-^cm\H,fi-^cm\E) eV{Mj^) = H(curr,fi) x H(curl,fi). 
Remark 3 

(i) Theorem \E holds if we replace the spaces by spaces of vector fields having such 
regularity only piecewise, i.e., C^, where $ G Cp, if and only if ^ e C^~^ and $ is 
piecewise C^. 

(ii) To obtain the second order regularity in Theorem\^ (iii) and in view of numerical 
applications it is sufficient to assume that {Eq,Hq) has H'^{Q)- components and that 
{F,G) has Q?{VLt)- components with bounded derivatives. 

If [E. H) admits the second order regularity of Theorem [2] (iii) then we can apply 
dt + iM^ to dSSI) and obtain 

{dl+Ml){E,H) = {F,G), 
where {F,G) := ((9* + i Ma)(F, G). Equivalently, we have 

{dtAdt-MA-^M){E,H) = A{F,G). 

Since 



curl /i ^ curl 

curl curl 
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the latter equation decouples for the electric field E and magnetic field H. 

In this paper, we intend to discuss the second order system for the electric field E, 
which reads in classical terms 

{dt e dt + cm\ fi'^ cm\)E = K := sF in ^t, (2.8) 

ux E\r = on Tt, (2.9) 

E{0) = Eo in n, (2.10) 

dtE{0) = E'o:= e'^cnil Ho + F{0) in fi. (2.11) 

Moreover, we assume throughout this paper that the second order regularity of Theorem 
M (iii) holds. 

Remark 4 A solution of the second order problem f l2.8p -f l2TTT|) provides also a solution 
of the original first order problem f l2.ip -f l23|) . In particular, under proper regularity as- 
sumptions on the data the system fl2.8p -f l2rTT]) is uniquely solvable as well. To show this 
it sufficies to set 



H(t) := [ {G{s) - fi-^cm\E{s))ds + Ho. 
Jo 



Then, f l2.3p and f l2.4p hold and fl2.2p and (12. 5p follow directly. Furthermore, to prove 
f l2.ip we use f l2.8p and the above definition of H and obtain 

dtE{t)= [ dl E{s)ds + E'q = -e-^ cml [ fi'^ curl E{s)ds + f F{s)ds + E'^ 
Jo Jo Jo 

= e-^ curl H{t)+ [ {F{s) - e'^ curl G{s))ds + E' - cml Ho . 

Jo ' V ' ^ V ' 

=dsF{s) =F{0) 

Hence, our further analysis is based on fl2.8p - fl2.lip . 

3 First form of the deviation majorant 

Let E be an approximation of ii^. In this section, we assume that 

C;(7,H(curl,fi))nC;(7,H(fi)). (3.1) 
Our goal is to find a computable upper bound of the error 

e:= E- E 

associated with E. For all t G / and p G we define two nonnegative quantities 

n^At) ■■= \ml,-^,p,nit) ■■= \\9t nln (t) + Pl^^^mU, nit) ^ (3.2) 
N^At) ■■= \ml,-^,P,n. ■= nln. + P llcurl , (3.3) 
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which generate natural energy norms for the accuracy evaluation. We note that by Fubini's 
theorem J^^ = and thus 



nq,^p{s)ds. 



Theorem 5 Let p G (0,1) and E he an approximation satisfying (13.11) . Moreover, let 
9t e G H(curr, Vt) for all t & I. Then, for all t ^ I 



where 



(3.4) 



-■ys r_ 



B 



E,p;Y,-y 



it) 



-7s 



E,Y,-t,P 



s)ds 



and the infima are taken over 7 G M+ and Y G Cp(/, H(curl, Q)). Here, 



fE,Y,-y,p 



9E,Y,-y,p ^-J 



,f,p 



7 



-1 



7,p ' ~E,Y^ 
2 



K 



E,Y 



dtK 



t ^^E,Y 



+ (1-P) 



-1 



K 



EX 



it), 



E,Y 



ne,i(0) + 2(K^5.,curle^^(0), 



where 



K 



EX 



K 



a?£,curiy 



ed'^, E + c\xi\Y - K, 



^E,Y ■= ^^^,iE,Y ■= ^ curl E-Y. 

Remark 6 We outline that the functional fEY-yp' ^Ep Y-y '^'^'^ ^Ep Y-y depend only on 
known data, the approximation E, the free variable Y and the free parameters p, 7, and do 
not involve the unknown exact solution E. Thus, these quantities are explicitly computable 
once the approximate solution E has been constructed. We note that the 'zero term' z^y 
represents the error in the initial conditions. In particular, z^y "mainly consists of first 
order derivatives of the initial error e(0) = Eq — E{0). Furthermore, Y may even be 
chosen from the larger space 

cl(i,H{n))ni^{i,H{cnT\,n)). 

Remark 7 The absolute value of the zero term 



Ze,y = \\9t e\\l^ (0) + fcurl eg_i_f^ (0) + 2 , curl e^^ (0) 
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may be estimated from above by the quantities 

Ze,y ■= ^e,i(0) + 2| (^k^Y, cnile'^ (0)| 

= \\dte\\l^ (0) + ||curle||^_.^^ (0) + 2| curie )^ |(0) 



e\\l^ (0) + 2||curle||J_,_^ (0) 



/i,f2 



(0), 



which are nonnegative and easily computable. Furthermore, 

(i) 2;^ y = 0, if dt e(0) = and curl e(0) = 0. 

(ii) % y = 0, if any only if dt e(0) = and curl e(0) = 0. 

(iii) %y = 0, if any only if dtcifS) = and curle(O) = and fiY{0) = curli?(0). 

(iv) QEx-y^p = 0, if and only if fiY = curl E and curl Y = K — e E. 

Therefore, choosing the functional fEVy p '"'^^^ ^ey ^"^ ^ey ^^^^ /^'^ all t E I the 

functional b^^.y^it) vanishes, if and only if 

dtE{0) = E'o, curl ^(0) = curl Eo, ^r = curlE, curlY = K - e E. (3.5) 

Thus, dfC and curie vanish, if and only if n^^p = 0, which is implied by ^.y^ = 0. The 
latter condition is equivalent to (13.51) . The same holds for the energy norm Ng^p and the 
functional Be p-y.^ 

Proof of Theorem [5] We start with deriving first order ordinary differential inequalities, 
which then lead to the estimates by Gronwall's lemma (see appendix). Since dt e belongs 
to H(curr,fi), we have 

dtne,i{t) = 2 {edfe,dte)^ (t) + 2 (/^"^ curl e, curl 5* e)^ (t) 

= 2(k -edfE,dte) (t) - 2 (^i-\m\E -Y + Y,cm\dte) (t) 



= -2 {^K^ y, dt ej^ (t) - 2 dt ^K^ y, curl ej^ (t) + 2 {^dt K^ y, cmlej^ (t). 
Thus, by integration 

ne,i(t) = z^y - 2( (K^y, curie) (t) - ((9t y, curie) + iKEy.dte) ). (3.6) 



=.Si 



--■■S2(t) 



■■S3{t) 



We estimate the scalar products Se by 

2|^i(t)| < a ||curle||^_i^^ (t) + 

2|52(t)|</?||curle|P... 



Kex 



dt Ke^ 



2 

p,Vtt 



it), 



(3.7) 



2\S,it)\<j\\dte\\l^^ + r' 



K 



EY 
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where a, /3, 7 G R+ are arbitrary. For p G (0, 1) we choose a := 1 — p G (0, 1). Then, for 
arbiratry 7 G R+ we put (3 := jp E Inserting fl3.7p into (13.61) . we achieve 

ne,p < lN^,p + fE,Y,-y,P^ 
which may be written in two ways 



Jo 

Gronwall's inequahties, i.e., (IA.2P and flA.4p . complete the proof. □ 



Since (^^ e = imphes e{t) = e(0) = Eq — E{0) constantly for all t we obtain: 

Theorem 8 Let an approximation E and Y as in Theorem [3] be given. Then, the 
following two statements are equivalent: 

(i) E{0) = Eo and b^^^.y.^ = 0. 

(ii) E = E and pY = curli?. 

In words: Let the approximation E satisfy the first initial condition E{0) = Eq exactly. 
Then, the functional b^ ^.y^^ vanishes if and only if the approximation E equals E and pY 
equals curlii^. 

Remark 9 The assertions of the latter theorem remain valid if we replace b^ ^.y^ by 

Remark 10 The latter theorem provides a new variational formulation for the second 
order problem fl2.8p - fl2.lll) and thus, in view of Remark\^ for the original first order 
problem fl2.ip - fl2.5p as well. 

3.1 Refinement of the estimate 

We can derive sharper estimates if p and 7 in Theorem O depend on time. Then, we 
replace 



ft 

n^,p{t), N^,p{t) = I n^As)ds 



^0 



[\\\dtm^{s)+p\\cm\nl^,^^is))ds 
Jo 



by 



:= Wdtnl^ {t) + p{t) llcurl*!^.,^^ (t). 



t ft 

2 



N^,pnit)-= / l{s)n^,p{s)ds = / As){\\dtm^^{s) + p{s)\\cni\^\^-,^^{s))ds, 
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respectively. In this case, N'^ ^^^ = ■yrii^^p and we modify f l3.7p in an obvious manner, i.e.. 



2|5i(t)| <«(t) ||curle£_.o (t) + a-\t) 



Ke,y 



2\S2it)\< / /3(.)||curler_.o(.)rf.+ / r\s) 



dtK 



E.Y 



{s)ds, (3.8) 











Jo 



{s)ds. 



By dniSD and ([3SD we find that 
Jo 

where fE,Y,'y,p '■= 9E,Y,'y,P + ^e,y with 



Ke,y 









Ke,y 


Jo 



2 



it) 



is)ds+ / i^p)-\s) 



dt K^ y 



{s)ds. 



We apply flA.ip and flA.31) . respectively, and arrive at the following result: 



Theorem 11 Let p : / — )■ (0, 1) and E he an approximation satisfying (13. ip . Moreover, 
let e G H (curr, d) for alltel. Then, for all t el 



ne,p{t) < ^B,p;y,7W' ^e,p,^{t) < inf B^^p.y,^{t), 



(3.9) 



where 



-r(s) 



7(s)/B,y,7,p(s)(is + /^y_^_p(t), T{t):= -f{s)ds 



e '^'^7(s)/s,y,7,p(s)(is 



.r(i) / «-r(s) 

'0 

anc? t/ie infima are taken over 7 : / — t- M_|_ and y G Cp(/, H(curl, Q)). 
Remark 12 

(i) The corresponding other assertions hold like in Theorem\^ 
(ii) 7/7 is constant we get the same formulas as in Theorem\^ i.e.. 



BE,p;Y,,it) = 7e^' / e-''^fE,Y,,Mds. 
Jo 

We note that in this case A^e,p,7 = 7A'e,p,i- Ifl and p are both constant the estimates 
coincide with those of Theorem\^ 
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(iii) Since the latter estimates are stronger it is clear that Theorem\^ and Rennarks\^ and 
[23 hold as well. 

4 Second form of the deviation majorant 

The estimates presented in Theorems [5] and [TT] are derived for approximations E 
having second order time derivatives. This requirement may be difficult to satisfy in 
practice because typical approximate solutions possess only ffist order time derivatives. 
In this section, we derive estimates applicable for approximations of such a type. 

As above, E is an approximation of E, but now we also introduce a vector field Et 
considered as an approximation of dt E. Hence, we define both the error and the error of 
the time derivative separately by 

e:=E-E, et:=dtE-Et. 

We note that in general Et ^ dt E and therefore Ct 7^ dt e. Henceforth, we assume that 

^,^,GC;(J,H(curl,fi)), 

et e H(curl^^]) for alH G /, 

where it would be sufficient to assume that 

E eCl(l,H{cuT\,n)), 
Et e Cl(I,HiQ)) nl\l,H{cm\,n)), 
Et e H(curr,r]) for all t G /. 

With two nonnegative, real functions p and 7 on / we define two energy norms 

n<p,^At) ■■= Wnln it)+p{t) ||curl^^||J_,^^ (t), 

iV$,*,P,7W := / 7(s)n<i,,vi,,p(s)ds = / 7(5) ( I'l'lJ (s) + p(s) 
Jo Jo 

Then, A^4,*,p,7 = T'^*,*,^- 

Theorem 13 Let p : / — (0, 1) and E be an approximation satisfying (14. ip . Then, for 
all t el 

ne,,e,p{t) < inf h,EuP;Y,'r(^)^ N,^,e,p,yit) < inf ^E.^^pjy.^l^). (4-2) 

where 

h,EuP;Y,,(t) e'^*^ [^-'^^'h{s)fE,E.,y,,,p{s)ds + fE,E,,y,,,,{t), r(t) := [\{s)ds, 



BE,EuP;yn(t) ■■= e'^*^ / ^~^^'h{s)fE,E.,y,,,,{s)ds 
Jo 
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and the infima are taken over 7 : J — )■ and Y G Cp{I, H(curl, Q)). Here, 



Key 



{t) + 2{K^y,cm\{Et-dtE) 



+ / 1' is)\\K^^y\\ {s)dS 

Jo 

+ I {lP)~\s) fi~'cm\ Et-dtY ' {s)ds, 
Jo M.f^ 



^E,Et,Y ■= '^et.e.ilO) + 2(^K^_y,curle^^(0), 
where K^^ y := /^a,^t,curiy ■= £ dt Et + curl y - K. 

Remark 14 // Ef = dtE then the estimates coincide with those of Theorem [711 Fur- 
thermore, Remark\^ holds in a similar way. Particularly, Y may he chosen from the larger 
space 

C;(7,H(^]))^L2(/,H(curl,^])). 

If > is constant then 

Net,e,p,jit) = -fNe,,e,p,l{t) = 7 / ^et,e,p(s)rfs, T{t) = -ft 

Jo 

and the upper bounds simply fy, i.e., 

bE,Et,p;Y,yit) = 1^^' [ ^~""fE,EuY,^,p{^)ds + fE,Et,Y,^,p{^)^ 

Jo 

BE,E,p;Y,,it) = 7e^* f ^-''fE,EuY,,^s)ds. 

Jo 

If both 7 and p are constant we have 



In this case, 



and 



N^,^,,{t) = \ml^^+p\\cnT\nl.,^^^ 

n<i.,*,p = iV$,^.p, N^,'f,pit) = / n.s>,^,pi^)ds. 

Jo 



Deviations from Exact Solutions of Maxwell's Equations 



13 



Theorem 15 Let p G (0, 1) and E he an approximation satisfying (14. ip . Then, for all 
tel 



where 



1 1 



Jo 



(4.3) 



^E,Et,p;Y,j 



it) := e^* / e-^7^,^,y,,,p(s)c^5 





and the infima are taken over 7 G M+ and Y G Cp(/, H(curl, Q)). Here, 

fE,Et,Y,l,P 9E,Et,Y,l,P + ^E,Et,Y-> 



9E,EuY,'r,p(^) ■= (1 - P) 



-1 



it)+2{kj,y,cnT\iEt-dt E) 



p,^t 



Remark 16 If Et = dfE then the estimates coincide with those of Theorem\^ Again, 
Remark Idi holds in a similar way. In particular, Y can be chosen from 

Cj(7,H(n))nL2(/,H(curl,fi)). 

Remark 17 The absolute value of the zero term 

^E,Et,Y = ll^tlle.n 

(0) + llcurl e||^_,^^ (0) + 2 (^K^ y, curl e)^ (0) 
can be estimated from above by the two quantities 



E,Et,Y ■~ '^et 



,1(0) + 2| ( K^y, cnile 



edl. o(O) + l|curle|M..o (0) + 2| K.^, curie) |(0) 



< z 



E,Et,Y ■- 



WetWln (0) + 2 llcurl ef„_.o (0) 



K 



EY 



(0), 



which are nonnegative and easily computable. The same manipulation can be done with 

the term 2 (K^ y, curl(£'t — dtE)\ , taking, e.g., it's absolute value, which leads to some 

\ ' 'fit 
nonnegative gE,Et,Y,'y,p^ I^oreover, 

(i) ^^y = 0, if et(0) = and curl e(0) = 0. 

(ii) % y = 0, if any only if et{0) = and curle(O) = 0. 
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(iii) -SgE^y = 0, if any only if et{0) = and curle(O) = and fiY{0) = curl£'(0). 

(iv) (jj^ EtY-y p~ 0' ^/ ^'^^ ^''^^y ^/ ~ "^^^l ^> ~ "^^^l ^'^^ "^^^l y = K — edtEt. 

Therefore, choosing the functional fE^EtX-r^p ^^^^ ^e,y %,y '^'^'^ 9E,EtX,'y,p ^^'^^ 
/or all t & I the functional bE,Et,p;Y,'yi^) vanishes, if and only if 

EtiO) = E'Q, curl ^(0) = curl Eo, (4.4) 

fiY = curl E, dt iiY = curl Et, curl Y = K - e dt Et. (4.5) 

Thus, et and curie vanish, if and only if net,e,p = 0, which is implied by bj^ Etp-Y-y ~ ^■ 
The latter constraint is equivalent to 04.41] and 04.51] . The same holds true for the energy 
norms N^,,e,p, ^et,e,p,7 and the functionals B^^^^^^.Y,-y' B^^Et,p;Y,i- 

Proof of Theorem 1131 We follow in close lines the proofs of Theorems 151 and fTTl Since 
et G H (curl", VL) and dte = et + Et — dt E, we have 

dt net,e,iit) = 2 (e (9t e*, et)f^ (t) + 2 curl e-Y + Y, curl Ct)^ (t) 
+ 2 (^fi-^ curl e, curl(^t - dt E)^^ (t) 

= 2(^K- curl Y ~edtEt,et'^Jt)-2 (^K^ y, curl e*^^ (t) 

+ 2 (cmliEt - dt E),fi~^ curl e^^ (t) 
= -2 {k^^^y, et)^ it) - 2 dt (k^^y, curl e)^ (t) 

+ 2(kj^y,cm\{Et-dtE)) (t) 



n 

+ 2 { /i"^ curl(Ef -dtE) + dt k^ y, curie ) (f). 



Thus, by integration 

"et,e,l(t) 



^^^^^y + 2 (K^^y, cm\{Et - a* i^;) )^ (4.6) 



2((K^_^, curie) (t) - ^/i"^ curl 9* F, curl e\ + (K^^ et)^J . 



= :S3(i) 



= :^l =:S2(t) 

If Et = dtE then 04.61) coincides with 03. 6p . As before, we choose a := 1 — p and /3 := 7p 
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and estimate the scalar products Se as follows: 



2\S,{t)\ < a{t) icurlei' o (t) + a-\t) 



K 



E,Y 



it) 



2\S2{t)\< / P{s)\\cnT\ef ,^{s)ds+ / ^'\s) cm\ Et - OtY {s)ds (4.7) 



2\S^{t)\< [\{s) WetWl^ {s)ds+ [\-\s)\\K^^^yf is)ds 
Jo Jo 

Inserting (14.71) into (14. 6 p yields 



ne,,e,p(t) < / 7(s)ne,,e,p(s)(is + /^_^^ y_^_^(t), 

Jo 

Finally, Gronwall's inequalities, i.e., (lA.ip and (lA.Sp . prove the assertions. 



Since Cj = implies 



we obtain: 



e(t)= / iEt-dtE)is)ds + e{0) 
Jo 



□ 



Theorem 18 Let approximations E, Et and Y as in Theorem [73 or Theorem [73 be 
given. Then, the following two statements are equivalent: 

(i) £^(0) = Eq and Et = dtE and h^^^^ ^.y^^ = 0. 

(ii) E = E and Et = dtE and jjY = curl E. 

In words: Let the approximations E, Et satisfy Et = dtE and the first initial condition 
£'(0) = Eq exactly. Then, the functional ^.y^ vanishes if and only if the approxima- 
tion E equals E and fiY equals curli?. 

Remark 19 The assertions of the latter theorem remain valid if we replace h^^^^ p.y .^ by 

Remark 20 The latter theorems provide new variational formulations for the second 
order problem (I2.8p - (I2.11I) and thus, in view of Remark\^ for the original first order 
problem (I2.ip - (l2.5p as well. 
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5 Estimates for the approximation of the whole so- 
lution 

By (12.61) (or the basic equations (12.11) . (I2.2p ) we also get estimates for the errors h, 
hf of the magnetic fields H, dtH and their approximations H, Ht- E.g., by adding 
-{Et, Ht) + iMA{E, H) to dlSD we obtain 

{ct, ht) - i MA(e, h) = if, g) := {F, G) - {Et, Ht) + i Ma{E, H), 

which reads explicitly 

et — curl h = f = F — Et + curl H, 
ht + curl e = g = G — Ht — /^^^ curl E. 
Therefore, we can estimate 

UhtApit) ■= pit) \\ht\l^n (0 + l|curl/i||J-i_f^ it) 

<2ne,e,p(t)+2||/||;_^ W + 2|bll',n {t) 

< 2mf + 2 \\f\l^ (t) + 2 \\g\\l^ (t), 

which yields 

net,e,p{t) + nh,,h,p{t) 

it) + pit) \Min it) + \\c^^^h\\e-^,n it) 
< 3mf + 2 (t) + 2 \\g\\l^ (t). 

Of course, similar estimates hold for the other norms and functionals and the estimates 
simplify in an obvious way if p or 7 are positive constants. The vector fields (/, g) measure 
the error in the original first order equation (12. 6p . Moreover, (e^, ht) may be replaced by 
dt{e, h) if for the approximations sufficient regularity is available. In this case, the error 
in the first order equation is 

{f,g) = (F, G) - {dt-iMA){E,H) = {dt-iMA){e,h). 

A Appendix: Gronwall inequalities 

Gronwall inequalities can be found in almost any book about ordinary differential 
equations. Since these estimates differ in small details, we present here two estimates, 
which meet our needs. With 

C;(7):=C;(7,M), C°(7):=C°(7,M) 

we have 
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Lemma 21 (differential form) Let u e Cp(/) and G C°(/) with > 0. If in I 

u' < (pu + ip 

then for all t & I 

u{t) <ex-p{^t))(u{0)+ [ exp(-$(s))V^(s)ds), := [ if{s)ds. (A.l) 

Jo Jo 

If If is a nonnegative constant then for all t & I 

u{t) < exp{(ft){u{0) + [ exp{-(ps)ijj{s)ds). (A.2) 

Jo 

// ^ < c e M then for all t el 

u{t) < {u{0) + ct) exp($(t)). 
Proof Since 

(exp(-<l>)M)' = exp(-$)(-u' - ifu) < exp{-^)'ijj 
we have ^ 

exp(-$(t))u(t) < u(0) + / exp(-$(s))^(s)rfs, 

Jo 

which proves the first part. The other assertions are trivial. □ 

Using the notations of the latter lemma, we have 
Lemma 22 (integral form) Let u,ip,'^ e C°(/) with <^ > 0. If for all t e I 

u{t) < [ ip{s)u{s)ds + ipit) 
Jo 

then for all t e I 

u{t) < exp^t)) [ exp{-^s))^{s)ip{s)ds + ip{t). (A.3) 
Jo 

If is a nonnegative constant then for all t & I 

u{t) < (f exp{ipt) / exp{—ips)ip{s)ds + ^|J{t). (A.4) 
Jo 

Ifip<ceR then for all t el 

u{t) < cexp($(t)). 
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Proof Set u{t) := / (p{s)u{s)ds and tp := (ftp. Then u' = (fu < ipu + ip. Lemma [5T] 
yields 

< + V'(t) < exp($(t))(M(0) + /" exp(-$(s))V^(s)rfs) + V'(t), 

which proves the first assertion since -11(0) = 0. The second assertion is triviaL To prove 
the last part we compute with exp(— $)' = — exp(— 

-u(t) < cexp($(t)) / exp(-$(s))v2(s)cis + c = cexp($(t)). 
Jo 



Remark 23 The differential and integral form are equivalent. 
Proof Let the assumptions of Lemma be satisfied. Then, 

u{t) < j ip{s)u{s)ds+ [ i!{t)ds + u{0) 
Jo Jo 



and by Lemma [22] 

< exp($(t)) / exp{-^{s))ip{s) ^{s) ds + ipit) 

Jo ' V '^-v^ 

t 4 

/•* . ^ t ^ 

= exp($(t))(/ exp(-$(s))^'(s)ds-exp(-<l>(s))^/'(s) )+tp{t), 







■ 



=exp(-${t))V(t)-V(0) 



□ 



which completes the proof since ^/'(O) = u{0). □ 



Corollary 24 

(i) Let the assumptions of Lemma l21\ by satisfied. Then, u{0) < and ip < imply 
u < 0. Hence, if u > then u{0) < and ip < imply u = 0. 

(ii) Let the assumptions of Lemma l2E by satisfied. Then, ip < implies u < 0. Hence, 
if u> then ip < implies u = 0. 
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